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ABSTRACT 

It is well known that our motion with respect to the cosmic microwave background (CMB) 
rest frame introduces a large dipolar CMB anisotropy, with an amplitude cc /? = v/c ~ 10~ 3 . In 
addition it should lead to a small breaking of statistical isotropy which becomes most notable 
at higher multipoles. In principle this could be used to determine our velocity with respect to 
the CMB rest frame using high angular resolution data from Planck, without directly relying 
on the amplitude and direction of the CMB dipole, allowing us to constrain cosmological 
models in which the cosmic dipole arises partly from large-scale isocurvature perturbations 
instead of being fully motion-induced. Here we derive simple recursion relations that allow 
precise computation of the motion-induced coupling between different spherical harmonic 
coefficients. Although the lowest order approximations for the coupling kernel can be deficient 
by factors of 2 - 5 at multipoles / ~ 1000 - 3000, using our results for the aberration kernel we 
explicitly confirm that for a statistical detection of the aberration effect only first order terms 
in p matter. However, the expressions given here are not restricted to /3 ~ 10~ 3 , but can be 
used at much higher velocities. We demonstrate the robustness of these formulae, illustrating 
the dependence of the kernel on [}, as well as the spherical harmonic indices / and m. 
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1 INTRODUCTION 

The large dipolar temperature anisotropy of the cosmic microwave 
background (CMB) is usually interpreted as a consequence of our 
motion with respect to the CMB rest frame, implying that the Solar 
System is moving with a velocity of f3 = v/c = 1.241 x 1(T 3 in the 
direct i on (l,b) = (264.14° ± 0.15°, 48.26° ± .15°) dSmoot et all 
Il977l : iFixsen et alj|l996l : IScott & Smootll2010l) . However, in ad- 
dition our motion should lead to a small breaking of statistical 
isotropy, which is most notable at high multipoles / or equiva- 
lently small angular scales, leading to the coupling of neighbouring 
spherical harmonic coefficients as a consequence of t he aberration 
and boosting effect dChallinor & van L eeuwen 2002j). In the direc- 
tion of our motion, the temperature anisotropies are beamed to- 
wards each other, while in the opposite direction they are magnified 
(see Fig.|2]for illustrati on). This should lead to a ti ny power asym- 
metry on the CMB sky dBurles & Rappaportll2006l) , and introduces 
correlations between neighbouri ng spherical harmonic coefficients . 

It was recently argued dKosowskv & Kahniashvilil l20ld : 
lAmendola et"a!]|2010h that the latter effect could be used to de- 
termine our velocity vector with respect to the CMB rest frame 
using high angular resolution data from the Planck satellite, with- 
out directly relying on the amplitude and direction of the CMB 
dipole. This is because the CMB provides both the most cosmo- 
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logically distant and the most statistically isotropic 'marker' on the 
sky, and hence can be used to search for the small aberration effect. 
In principle this should allow us to constrain cosmological models 
in which the cosmic dipole a rises partly from lar ge-scale isocur- 
vature perturbations (e.g., see IZibin & Scot3l2008l and references 
therein) instead of being fully motion-induced. Alternatively, one 
could use the small motion- i nduced asymmetry in the S Z cluster 
dZeldovich & Sunvaev|[T969l : ISunvaev & Zeldovichlll980l) number 
counts to independently measure our velocity with respect to the 
CMB and place constraints on the primordial dipole. However, the 
required number o f SZ clusters is rath er large, which makes this 
endeavour difficult dChluba et al.l2005l) . Also in the future, the cos- 
mological 21cm signal from the reionization epoch could provide 
another opportunity to search for the aberration effect. 

To account for the effect of our motion on the CMB 
anisotropies one has to compute the amount of mixing between 
neighbouring spherical harmonic coefficient^ a lm . Here we restrict 
ourselves to the CMB temperature anisotropies, however, it should 
be possible to extend our method to the case of polarization. It was 
shown earlier that in lowest order of (3 our motion leads to a cou- 
plin£0of «;„, with a/±i„>. However, it is difficult to compute the aber- 



1 Alternatively, one can directly work in real space, however, here we 
follow the example of earlier works ([ Challino r & van Leeuwenl |2002| ; 
iKosowskv & Kahn iashvili 201cl; lAmendola et all20icl7 on this problem. 

2 This assumes that the z-axis is aligned with the velocity vector. 
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ration kernel which describes this coupling, since the associated 
integrals are highly oscillatory, making numerical quadrature very 
demanding and time-consuming, even for the lowest order coupling 
terms, i.e. between / «-» / ± 1. One way around this problem is to 
use series expansions of the corresponding integrals in orders of 
P <k 1. Analytic expression accurate up to 0{f} 2 ) were obtained 
earlier, h owever, at large / these expressio ns converge slowly once 
ip > 1 dChallinor & van Le euwen 2002). Furthermore, in recent 
discussions of the aberration effect only the lowest order expres- 
sions, i.e. 0(fS), were applied. As we show here, these can be defi- 
cient by factors of 2-5 at / ~ 1000- 3000, and it becomes important 
to include higher order terms when computing the coupling kernel. 
However, for the statistical properties of the CMB only the first or- 
der terms in p really matter, as we explicitly confirm here using our 
results for the kernel (see Sect. [3). 

For this purpose, we derive general recursion relations that 
allow accounting for terms up to high orders in p. Rather than fo- 
cusing on p ~ 10~ 3 , we discuss the kernel for more general cases, 
showing that our method is both very fast and very robust, even for 
much higher velocities. We find that for precise computation of the 
aberration kernel at multipoles / > 1000 - 3000 in the case of our 
motion with respect to the CMB rest frame terms up to high orders 
in p (e.g., O0'~'' l p s ) at / ~ 3000) are required. We also illustrate 
the I, in and p dependence of the aberration kernel. 



2 THE EFFECT OF (OUR) MOTION ON THE CMB SKY 

In this section we briefly recap the key formulae to take the ef- 
fect of motion on the CMB temperature anisotropics into ac- 
count. Our formulation most c losely reassembles the one in 
iKosowskv & Kahniashvil] feoicl ). however, here we do not use a 
first order series expansion of the problem, but give general recur- 
sion formulae that in principle allow us to compute the coupling of 
different modes to machine precision for practically any value of p 
(see Appendix lAl. 



where Q(x) = xe"/[e x - 1] and xi = hv'/kT^. Here T' Q defines 
the apparent monopole temperaturqfj in S' , and V denotes the fre- 
quency at which the measurement is performed. In addition, we 
have assumed that A' « 1. However higher order corrections could 
in principle be taken into account, in next order leading to a y-type 
spectral distortion when compari ng with a blackbody of tempera- 
ture T' l lChluba & Sunvaevll2004l) . 

Rotating all z-axes parallel to the velocity vector of the moving 
frame one simply has <p = <f>' and ji = +/3]/[l + /?//], where 
yti = cos(0) and p! = cos(0'). The expression for T and A' then read 

r'(0',f) = r u [i + AX0',f)], 

T 1 +A(0,0') 



- 1, 



(3a) 
(3b) 



where we used v/T = v'y[l +Pfi']/T = V IT' , In addition, should 
be expressed as functions of 0' using the relations from above. To 
simplify matters further we write 

AV, 0') = A' (fi', (/>') + A^Qi', 0'), (4a) 



1 



To A(ji,(f>) 



- 1, 



(4b) 



(4c) 



T'o yD +m 

Here A' (ji' is the temperature anisotropy in the moving frame 
that is arising from the CMB rest frame monopole term, oc T , 
alone. Due to Lorentz-boosting it results in a motion-induced 
dipole, quadmpole, octupole, and higher order multipoles, all with 
increasing order of p. Since /3 <sc 1, for our Solar System one usu- 
ally can stop after the motion-induced quadmpole. The contribution 
A^ n (/j', (p') arises from the primordial CMB anisotropies. Here both 
boosting and aberration terms are contributing, with the aberration 
terms dominating at small scales. 

One can now perform a spherical harmonic expansion of the 
temperature field, T'(6',cf>'), to obtain the spherical harmonic co- 
efficients a' ln: (P) that describe the CMB sky inside <S'. This yields 



2.1 Basic formulae 

In the CMB rest frame, <S, the energy spectrum of the CMB black- 
body is given by /,,(#, (p) = B V (T), where B v is the blackbody spec- 
trum of a given temperature at frequency v, and T = T(6, <p) de- 
scribes the CMB temperature in different directions of the sky. Fur- 
thermore, one can write 



T(6>,0) = r [l+A(M)], 



A(M) = ^^<i/»W^), 

;=i m=-( 



(la) 



(lb) 



using a spherical harmonic expansion of the temperature field. Here 
To is the value of the CMB monopole in S and A describes the pri- 
mordial CMB temperature anisotropies, including a possible dipole 
asymmetry created, e.g., by large-scale isocurvature perturbations. 

To transform 7 r (0, (p) into the moving fram^l S' we use the 
Lorentz-invariance of the photon occupation number to obtain 



i'A0',<P') 



2hv' : 



1 



c 2 ghv/kT _ J 



*fl v X^)[l+0(y o )AX0', </>')], (2a) 



3 At high / the aberration effect dominates over the Doppler term. However, 
the aberration kernel as defined here includes the contributions from th e 
Doppler term, which for the CMB statistics matters jAmendola et alj2010h . 

4 Henceforth primed quantities denote variables in the moving frame. 



AX0', ¥) = YjYj Ylm(ff > lp ' ) ~ 1 
4„(P) = J K ; >',0')AV^W 

CO /' CO 

= Z Z <;"'v3)^,„^X<" ,( ^ )a ''' 



(5a) 



(5b) 



The aberration kernel %\"' (P) is defined in Appendix lAl For the 
last equality, we have used the fact that in the chosen frame the ker- 
nel only mixes the different / and /' for fixed m (see Appendix lAl. 
The general case can be obtained by rotating the kernel. Note that 
«oo = V4^ by construction. 

Equation J5bt shows that the problem is fully determined once 
all elements of the aberration kernel are known. In Appendix lAl we 
explain how to obtain these using simple recursion relations. 



2.2 CMB monopole term 

To understand the effect of the motion of our Solar System on the 
CMB anisotropies we start by writing down the motion-induced 



Although for our Solar System it is clear that TL w Tq, one can derive 
more general expressions that account for the leakage of power from the 
higher multipoles into the monopole term. 
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Figure 1. Modulus of the aberration kernel < 7CS,°(jS) for different / and in- 



creasing value of /}. In the computation k 
reach convergence. 



256 was used for /? = 0.99 to 



signal arising from the CMB monopole alone. Although this case 
can be easily treated using a series expansion in terms of /3, it illus- 
trates the procedure in the more general case. 

To give the solution to the problem we need the kernel ele- 
ments "K^ifi). From Eq. JA14ab one directly has 



<>' 1=0 



p 



2k + 1 



roar^anh^To 
py t 



with ^k^ = l/[2fc +1]. This result can also be easily verified 
by analytic integration. To obtain the elements r K l ^(P) for / > we 
can use the recursion formula Eq. dA20t . With the initial condition 



1 it is straightforward to show that 



(0)-(M 
K 



(0)~0->/ 



(2/+ 1)!! 
g/(/ + 4)! 
8 (2/ + 5)!! 



(0)-0-K _ 



g/(/ + 2)! 
2(2/ + 3)!! 

g, (1 + 6)1 
48 (2/ + 7)! 



(7a) 



(7b) 



with g; = V2T+T. With these coefficients one can determine all 
"K^ifi) up to 12 th order in {S, while higher orders could be easily 
obtained with Eq. JA20b . Furthermore, according to Eq. (IA51 one 
directly has the kernel elements K'^(J3) = (-l)'K^(J5), which allow 
us to compute the leakage of power from the higher multipoles to 
the apparent CMB monopole. Putting everything together we find 



n 



arctanh(/3) 



6 



(8) 



P Qio 2/? 2 a 2 o 
VI ^An 3V5 V4tt 
for the value of the apparent CMB monopole in <S'. Clearly higher 
order terms are very small in the case of our Solar System, and as 
Fig.[T]shows, even for/3 ~ 0.9 only about half of the intrinsic dipole 
is expected to leak into the monopole term. Similarly, we find 



1+A'( 



!>')■ 



T arctanh(y6) 



j, 
1 (1 



py 



*1_ V^4= F io(^/)+ V^r-^-F^A (9) 
V3 3V5 

where the sum accounts for the motion-induced anisotropics aris- 
ing from the CMB monopole only (i.e. leakage of power from 
/' = to / > 0). Here we approximated T » T' Q . As is well 



known, to leading order only the motion-induced dipole really mat- 
ters in the case of our Solar System. However, the recursions given 
here also allow us to compute the effect for large /?. For exam- 
ple, a fast moving electron inside the hot gas of a galaxy clus- 
ter (with kT c ~ 5 - 15keV) can have velocities p ~ 0.1. Be- 
cause of the leakage of power from the monopole term to higher 
multipoles the electron 'sees' a reduced CMB monopole, while 
dipole, quadmpole and higher multipoles are increased accord- 
ingly. In the rest frame of the moving electron the scattering 
process is given by Thomson scattering, implying that only the 
value of the rest frame monopole and quadrupole matters. This 
change of the CMB monopole and quadrupole in the rest frame 
of the scattering electron is one reason for r e lativistic corrections 
dChallinor & Lasenbv||l998l; lltoh et al.lfl99il ISazonov & Sunvaevl 
120001: Chluba et alJl2005T) to the SZ signals from galaxy clusters 
jZeldovich & Sunvaevlll'969l ; ISunvaev & Zeldovichiri98c|) . 

Note that the change in the apparent monopole temperature 
also means that, as a matter of principle, one should correct the 
observed monopole to the CMB rest frame value when calculating 
the CMB power spectra. However, in practice this correction will 
be be low the cosmic variance level of the monopole JZibin & Scottl 
l2008h . and hence is negligible (D. Scott, private communication). 



2.3 Motion-induced terms for higher multipoles 

In this section we present several results for the aberration kernel, 
illustrating its dependence on /, /', m and /?. These were obtained 
using the recursion relations given in Appendix [A] however, we 
checked the precision of the results in several cases using explicit 
numerical integration, finding relative differences < 10" 10 when us- 
ing a large numbei^of terms in the recursions. 



2.3.1 Illustration of the aberration and boosting effect 

To illustrate the effect of motion on the CMB te mperature 
anisotropics, we modified Healpix dGorski et al.| |2005h to allow 
accounting for the aberration and boosting effects caused by the 
Lorentz transformation of the CMB sky into the moving frame. 
The results are shown in Fig. [2] We chose very high velocities here 
to clearly illustrate the effect. In the direction of the motion (up- 
wards) the anisotropics are beamed towards each other, while in 
the opposite direction the effect of aberration acts like a magnifying 
glass. One can clearly see a asymmetry between the power in the 
southern and northern hemispheres. Of course, for our Solar Sys- 
tem p ~ 10~ 3 and the effect is much smaller and can only be picked 
up by looking at the breaking of statistical isotropy. Also the power 
asymmetry will be much smaller. While the asymmetry is expected 
to be oc 2p when computing the difference between the power spec- 
tra on the two hemispheres, on t he full sky it is only of or der 1 
dChallinor & van Leeuwenll2002l juries & Rappaporul2006l) . 



2.3.2 Dependence on I' and p for in = 

In Fig. [3] we show the behaviour of the kernel, 7C ; ' (j6), with dif- 
ferent /' and p ~ 10~ 3 (upper panel) and for /' = 500 but varying p 
(lower panel). For clarity we have plotted the modulus of the kernel 
to suppress its alternating behaviour. As was expected, the kernel 



Depending on the value for f} this could mean k ■ 
k ~ 256 for large / and /}. 



4 for / < 3000. but also 
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Figure 2. Illustration of the effect of motion on the CMB temperature 
anisotropics for high velocities along the z-direction towards the north pole. 
The upper panel shows the CMB in its rest frame, the middle panel is for 
fl = 0.5, and the lower panel for /} = 0.9. In all cases the monopole and 
dipole terms were not included. The m aps were created using a modified 
version of Healpix iGorski et al.l2005h . 

becomes broader with increasing /' and fixed value of /?, and sim- 
ilarly, for fixed /' but increasing p. Also the modulus of the kernel 
appears to be quasi-symmetric, however, this symmetry is not per- 
fect, since there is a small but non- vanishing derivative of 'K'^ifj) 
with respect to /', which makes |*7f/,'° A , (j8)| t \'K' I ', AI0 (J3)\. For prac- 
tical purposes this quasi-symmetry could be used to compress the 
kernel and minimize storage. 

Figure [3]also shows that even for fi = 1(T 3 (i.e. close to the 
value of our own motion with respect to the CMB rest frame) at 
V > 1000 the coupling of F — > /' ± 2 becomes important. For 
/' ~ 3000 it is even stronger than the coupling /' — > /' ± 1, and 
also the /'—>/' ± 3 and /'—>/' ± 4 terms start to be significant. As 
mentioned above, this demonstrates that for /' > 1 //? higher order 
terms in the series become important. 

To push this point even further, in Fig. |4]we show the aber- 
ration kernel for ft = 10~ 2 . One can clearly see that the kernel be- 
comes very wide, rendering a perturbative analytic expansion diffi- 
cult. In the recursions we used k — 64, i.e. included terms £ = /'/? 
up to £ 128 . The computation for all elements, < 7C / ' °(/?), with V < I 
and |A/| < 100 up to /' = 3000, takes a few seconds using the re- 




-10 -8 -6 -4 -2 2 4 6 8 10 



A! 

Figure 3. Modulus of the aberration kernel 9fm°C8) for different Al = I - /'. 
Upper panel: for different /' and /? = 10~ 3 . Lower panel: for /' = 500 and 
different values of /J. 

cursion formulae, while brute force numerical integration just for 
/' = 3000 takes hours with Mathematica. We also implemented an 
integration scheme based on Chebyshev quadrature, however, also 
in this case the computation takes too long for real applications. 

2.3.3 Dependence on m 

In Fig. [5] we illustrate the m-dependence of the aberration kernel. 
For fixed /', the kernel becomes narrower with increasing m. Also 
it is clear that the variation of the kernel with m is rather slow. 
Apparently the important parameter is ml I. This fact allows a strong 
compression of the kernel functions, e.g., for /' = 2000 the kernel 
form ~ 0-500 is practically not changing. Therefore one only has 
to store the coefficients ''"'/c'^' for a fractions of the kernel elements, 
to obtain very precise results. Numerical computations can benefit 
from this property, allowing memory consumption to be reduced. 

2.3.4 Convergence of the recursions 

We already mentioned several times that for large /' and/or large f3 
many terms in the series have to be computed. This is because the 
coefficients i "' y K , -,' k * 1 strongly increase with k. In Fig. [6] we show the 
convergence of the aberration kernel for /3 = 10~ 3 at representative 
values of At /' ~ 750 the lowest order terms (i.e. 0{fi) for /' — > 
/' ± 1 ; O0 1 ) for /'—>/' ± 2, etc.) already give rather good results for 
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Figure 4. Modulus of the aberration kernel 9C£°(j8) for/J : 
V . In the computations we used k = 64 to reach convergence. 



10 2 and large 



750 



the strength of the mode coupling. We note that even for /' 
the kernel element for /'—»/' ± 2 has amplitude ~ 10%. 

Going to higher values of /' the kernel becomes broader, as 
explained above. For V ~ 2000 - 3000 terms up to O0'-'' 1 /3 s ) 
are important for accurate computations of the aberration kernel. 
The recursions given here easily allow such precision, while be- 
ing sufficiently simple. Furthermore, as our computations show, 
a lowest order expansion will not allow us to compute the pre- 
cise value of the coupling kernel once /' ^ 1000, as also 

mentioned earlier ( IChallinor & van Leeuwe i 120021) . We show 
this more quantitatively in Table [TJ where we compare the re- 
sults for the kernel obtained wi th the approximations given by 
iKosowskv & Kahniashvili t201Ch to our result for k = 6. The case 
o flKosowskv & Kahniashvilit201Ch is equivalent to our case k = 0, 
and we confirmed the values by direct comparison with their ex- 
pressions. As Table [T] clearly shows, the lowest order expansion 
can be deficient by a factor of 2 - 5. 



3 SIGNAL COVARIANCE MATRIX AND 
MEASUREMENTS OF /?. 

Although the values of the kernel elements computed to first order 
in p can be deficient by factors of 2 - 5 (see previous Section), for a 
statistical detection of the aberration effect it is important how the 
signal covariance matrix is affected. In real space the terms arising 
because of aberration can be modelled as dipolar convergence (e.g., 




Figure 5. Modulus of the aberration kernel l Kl' m (J3) for fixed j} = 10 3 and 
different /' and m. In the recursions we used k = 6. 



see lAmendola et al J2010h . In th is case, analogy to the CMB lensing 
effect ( ILewis & Challinoti r2006) implies that higher order terms in 
P will not be important for a statistical detection of the aberration 
effect, as order by order in leading terms canc el. 

As pointed out by lAmendola et alj d201fj|) at high / even the 
Doppler term matters. The Doppler term is not directly analogous 
to a lensing effect and hence the simple argument about its impor- 
tance for the signal covariance matrix is not evident. However, one 
can use 1/[1 + /?//] = YtZoi-^Wf in the definition of the ker- 
nel, Eq. JAU . Repeatedly applying the recursion relations for the 
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Figure 6. Modulus of the aberration kernel 9CL°(/S) for p = 1CT 3 and dif- 
ferent V . In this figure we illustrate the convergence to the series expansion, 
Eq. )A14at . For /' > 1000 several correction terms have to be taken into 
account to obtain precise results. 



Table 1. Representative kernel elements for /' — > /' and /'—>/' + 1. We 
compare the result obtained with the approximations Eq. (13) and (14) in 
iKosowskv & K ahniashvil1 [20 id) with those from this work up to k = 6. In 
some cases the lowest order expansions are deficient by a factor of a few. 
For /' > 2500 even the sign changes for the kernel element 70',^ (y3). 



1 




a /0 \p) in u\p) 


A /0 \p) in u\p ) 


1500 
1500 




1 


1 

-0.75050 


0.51155 
-0.55808 


2000 
2000 




1 


1 

-1.00050 


0.22360 
-0.57666 


2500 
2500 




1 


1 

-1.25050 


-0.04863 
-0.49685 


3000 
3000 




1 


1 

-1.50050 


-0.26021 
-0.33869 



rections in the c ovariance matrix, in agreement w ith earlier works 
on this problem (Challinor & van Leeuwen 2002). 

With our results for the kernel elements we were able to con- 
firm this statement using explicit computation of the covariance 
matrix elements (aj* aJ +lm V Although we explicitly use coupling 
terms up to A/' = 10 obtained with our recursions, the final covari- 
ance matrix element agrees with the simple first order result 



( a 'i>' Mm ) x -P'J + l ) 



a+D 2 



4(1+ l) 2 - 1 



[C M - Q] (10) 



where C; denotes the CMB temperature power spectrum. We found 
that all other off-diagonal covariance matrix elements with /' > / + 1 
for an ideal experiment do not affect the signal-to-noise ratio by 
more than ~ 0.5% - 1% for / < 3000, and hence can be neglected. 



4 CONCLUSIONS 

We have obtained simple recursion relations that allow computing 
the elements of the aberration kernel in a very fast and accurate 
way. The application of these recursions is not restricted to the 
small fi case, but in principle enable precise computations of the 
aberration effect for general p. Here we illustrate the main proper- 
ties of the aberration kernel for a wide range of parameters. 

Using the recursion relations we show that the lowest order 
expansions for the couplings between spherical harmonic coeffi- 
cients /' — > /' and /'—>/' ± 1 can be deficient by factors of 2 - 5. 
For /' ~ 2000 - 3000 terms up to 00''^ /3 s ) are important for ac- 
curate computations of the aberration kernel. The recursions given 
here easily allow achieving such precision, while being sufficiently 
simple. Albeit this large difference in the values of the kernel, for 
a statistical detection of the aberration effect the only the covari- 
ance matrix elements (fl£,aj + i m ) really matter, and the first order 
expansion in /? provides a sufficient estimate. 



associated Legendre polynomials allows absorbing each of the ex- 
tra factors of fx', order by order. In this way it is clear that also the 
Doppler term essentially lead to a coupling of neighbouring modes 
that is similar to an aberration effect, however, in each order of /? it 
is about / times smaller than the aberration effect. Still in the covari- 
ance matrix it contributes at the same level lAmendola et"aL 
Again by analogy with the CMB lensing effect this implies that 
also higher order terms from the Doppler effect lead to small cor- 
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APPENDIX A: DEFINITION AND COMPUTATION OF THE ABERRATION KERNEL 

The effect of our motion on the CMB temperature anisotropies can be fully characterized by the aberration kernel, which is defined as 

where \i = yu(//,/?) = (// + /3)/(l + /?//')• Using the definition for the spherical harmonic functions 

Y lm ( M ,<p) = N lm P';'(M)e""*, (A2a) 



Nl <" = \l — a 77T m ' (A2b ) 

y An (l + my. 

where -P™(/i) are the associated Legendre polynomials, it is easy to show that 

<T' 08) = ^ <W N lm Nv m ( A3a) 

^=L^rhir ldfl - (A3b) 

The problem is now that the explicit computation of the integrals r H",'^ l (f}) is a cumbersome task, since for large /, /' and \m\ <k /, /' one is 
dealing with highly oscillatory functions. However, as we will show below, it is possible to obtain simple recursion relations that avoid this 
problem and allow precise and fast computation of H",'^^) to machine precision for practically any value of p. 

Al Simple properties of the integrals "H/f^C/?) 

Knowing the properties of the associated Legendre polynomials, it is simple to deduce some useful properties of the integrals H",'^^), and 
hence the aberration kernel. From P"'(-x) = P"'(x) it directly follows that 

= (- 1 )' + '' 03). (A4) 

This implies that the de-aberration kernel r Kj n ' i " (-/}) can be directly obtained from % , h "' (/?). Furthermore, with the identity dyu'/[l + = 
dyu/[l one has 



= <Hl'l P = (- l) l+l '<H;i,, (J3). (A5) 

This relation reduces the number of independent kernel elements by a factor of two. For example, knowing how the CMB monopole leaks 
into the higher multipoles, i.e. knowing < H^ I (P), we conveniently have , H?_ >O (0) = (-l/'W/OS) for all /. 
Also, with Pj m (x) = (-l) m S P?W we obtain 

(I - mV (/' -mV 

nZ, 08) = jr-^- jj-^- 7C- 08), (A6) 
(/ + my. (/' + my. 

which implies that one only has to compute the aberration kernel for m > 0. Another compression is caused by the fact that the kernel 
elements are (steeply) decreasing functions with increasing Al = I - !'. Therefore in numerical applications, for every pair (l',ni) the kernel 
elements have to be computed until / = /' + A/ miix , since for Al > A/ milx the coupling elements vanish for practical purposes. For ji ~ 1CT 3 one 
expects A/ max ~ 4 - 5 to suffice at /' < 3000. 

Finally, according to Eq. dA3| l. the aberration kernel always acts on spherical harmonic coefficients with fixed m. It is also clear that 
< H",'^i{P) = for m > /', simply because there is no initial combination (l',m) — » (/, m) with in > /'. This property of the kernel simplifies 
the computations significantly. For example, the monopole can only leak into the m = terms of the higher multipoles, the dipole can only 
leak into the m = (-1,0, 1) terms of the higher multipoles, plus the m = term of the monopole, the quadrupole can only leak into the 
m = (-2, - 1, 0, 1, 2) terms of the higher multipoles, plus the m = (-1,0, 1) terms of the dipole, and the m = term of the monopole, and so 
on. 



A2 Recursions for the integrals r H",'^ l (J3) 

We now derive simple recursion formulae for the integrals < W ; "'^ ( (/3) for /' < /. It turns out that one also has to restrict the relations derived 
here to the case /' > in. The case /' = in will be discussed below (see Sect. lA2.2l >. We start with the recursion relation 

(/ + 1 - m) Pl.Qx') = (21 + Y)n' />>') - (I + m) Pli(M') (A7) 

for the associated Legendre polynomials. Adding (21 + l)/3 PT(jx') on both sides of the equation, and multiplying by PJ"(//)/(l + /?//), after 
integration over d// one obtains 
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The remaining integral can be simplified in the following way: 



- 1 -nw~ ** i + ^ d ^ ■ 

Now, using (21 + 1) xP'"(x) = (/ + 1 - m)P'" +l (x) + (I + m)P'" { (x) and rearranging terms, together with Eq. (IA8I we obtairQ 



(A9) 



(A10) 



Here we introduced the abbreviations a 1 " = [2/ - l]/[/ - m] and /!'" = a'" - 1. 

The problem is now that numerically the recursion Eq. (IA101 is not very stable, since leading order terms cancel each time the recursion 
is applied. However, one can use a Taylor series for Hf^, and derive recursions for the series coefficients instead. The parity of the kernel 
suggests the ansatz 



(All) 



Here we also used the fact, that the kernel coefficients are decaying with increasing M =\l- 1'\. Inserting this into Eq. l lAlOt . and rearranging 
terms, we obtain 



i-t/'-i _ 2(m) J /-i^r 



for / = /'. Here '"'V^' = for A: < 0. Similarly, for / = /' + 1 we have 



2a"r"'k 



(mU'-.;'+l 
K 2k 



In'" 



"T+2 ^ in I K 2*-2 K 2k-2 



A l' + l 

2a"' 



i» (my-i-»;'-i _ 7 m (my-i-.;' , /'+2 

'ty+l K 2Jfc z "/' + l K 2£ + „m K 2Jt-2 



?„'» [ ''+2 2*-2 2/ 



i-»/'+2| 

2A-4 



(m)„/' + l->/'+2 



2i-2 



Z, "/'+l a /'+2 



and for / > /' + 2 we find 



(my-./ _ _L h™ (")/-»/-! . W/-»/+il 

2 * ~~ /V" L 2* "■" K 2k-2 \ 

rr'" 

1 V \im (n!)/-l->/-l _ „i« (my-l-H , (my-l->/+ll 
, v «i [ '+I 2* "7+1 K 2t-2 ^ K 2*-2 

"/+1 

3 m 

+ -J—\]l m (i»y-2-»/-l _ „'» W/-2-»/ . (»)/-2-> 
m [ / +I 2 *-2 "7+1 "2/1-2 T K 2k-A 



(A 12a) 



(A 12b) 



(A12c) 



Equations JA12b in principle fully determine the problem and allow us to compute the aberration kernel in a fast way. However, two pieces 
are missing: (i) the initial condition *H™_, m ; and (ii) the relations for 'H™,. We derive these in the next two sections. 



A2.1 Initial condition for the recursions 

To compute the elements ( '"k^p m we use the identity P™(x) = (-\) m (2m - 1)!!(1 - x 2 )'" 12 and insert this into the definition of K"!^ v After 
some algebra, this leads to 

(Ai3) 

7 f^i 2" k\ (2m + 2/c+l)!! 

To define the initial conditions for the recursions, it is convenient to redefine the integrals — » y^""H"'^ r The recursions Eq ( IA12t 

remain completely unaltered by this transformation, but one can avoid performing another series expansion of y' m in p. Also, since the kernel 



7 For this equation we have assumed that /' > m, so that the case /' = m has to be treated separately. 
© 0000 RAS, MNRAS 000, 000-000 
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normalization factors are strong functions of /, /' and m (see Eq. JA2b ). it is useful to rescale the integrals *HT_^, in addition, absorbing the 
factor Ni'„, N/,,,. Afterwards, the aberration kernel can be expressed as 

<TOT = (-D' + '' P C£ f "* ) 4^ > (A14a) 

( '">4'f ' = 2nN Vm N lm (A14b) 

fM) (2k + m)\ (2m + 1)!! 

2t " 2*Jt!m! (2m + 2* +1)1!" ( } 

After this rescaling, the initial condition for k = is Wg!^"" = 1 . It is also clear that for in > 4 one always finds ( '" ) ^™^'" > 1 . In particular, 
for m » 1 and m » /; one has f m 'K~"" ss w*/[(4e)*fe!], such that for /? ~ 1CT 3 the series of c K"^ m (fS) always converges rather fast while 
m < 10 4 . Furthermore, it is directly cleai0 that (m lS^'~ > '' = 1 for all I', since for ft = the aberration kernel should be = 8m S,„> m . 

The rescaling also affects the recursion formulae Eq. (IA12I . One simply has to multiply all terms with appropriate ratios of 
[Nvm Nt m ]/[N im N jm ]. Defining a™ = ^/[4P - l]/[/ 2 - m 2 ] and 1™ = a^/a^_ v we find that for Eq. dA12t one simply has to replace a -> S 
and /!—>/!, without any further changes. 

A2.2 Recursions for fi 1 " , 

To obtain the required formulae for "H™, we start with the identity 



P^iW = -(2m + 1) VH?^W. (A15) 
Multiplying on both sides with />" !+1 (//)/[! + and integrating over d/j' we find 



)£p;- + V)^|^; 



-i 

where we made use of the identity dyu'/[l + Pi 1 '] = dyu/[l - Pfi]. Furthermore, 



-_ VT^7 



y[i-A»]' 

such that we have 

Xi ^ +V) ~i~^|~ p>)d/u = r X! ^'(p'lvn^^w^. (Ai?) 

Then with -y/l - // 2 Pf +1 (/z') = 0~ m)jj!P]*(j/) - (I + m) PJ_ x (ji') it is straightforward to show thaO 

y- 1 = (2m - 1) [(/ - 1 + m)^, +/3(/ + 1 - »)^] 

-/?(/- 1 + m) -TC" 1 ,. [ - (/ + 1 - m) «H£j . (A 1 8) 

Here again it is better to express this relation in terms of the series coefficients ^Tt^T* using the ansatz Eq. iJaTTJ. We again replace "H'^i = 

m—*l 



(m)~m->i = V2m +Ta [b e»-i)j^- w-i _ (™-D^-i 2 -/j 

+ a [/^'"-'Vjp'- 1 - fm - 1) ^'l + ( " ^-2. ( A19 ) 



where a = fe^l and = [|±± ^1 for m > and / > m. 

L 2m{l+m) J L 2/— 1 J 

In a similar way one can show that for m = 



<C) *?f ' = , l m St"- 1 + (A20) 



ViF^T V4/[/ + 2] + 3 

This equation closes the problem, and we are ready to compute all elements of the aberration kernel to machine precision. For this one can 
first generate all the required coefficients a'" and then compute order by order in k until convergence is reached. Looking at the properties of 
the recursion formulae Eq. JA12L JA19) and ( IA20K it seems easiest to use the following procedure: 

(i) Starting with '"fe^ for all k < k mdx one can first compute '"fe^' f° r a ^ required I > using Eq. l lA20t . 

(ii) Next determine ^°J^Z*' for all required / > 1, followed by < -°%%Z' 1 , '"fc 3 ^' until subsequently applying Eq. J A 1 2at — J A 1 2c b . 

(iii) Using Eq. l |A19t and the initial condition Eq. JA14ct compute all "fe-I^' for all required / > 0. 

(iv) Next determine {ly K^' for all required / > 2, followed by ^tk' untl l (1 fe2™ x ^' m ", subsequently applying Eq. l lA12at - jA12ct . 

(v) Repeat (iii) and (iv) for all required 2 < m < / max . 

8 This can be also easily shown by expanding 'H2 , _ ) ,, to lowest order in fi, i.e. < W7 t „ as J P"',(x) P"!(x) dx = 1. 

9 Here the implicit assumption is that m > 0. 
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